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Abstract. Following work of Tataru, 1151 and 1131 . we solve the division prob- 
lem for wave equations with generic quadratic non-linearities in high dimen- 
sions. Specifically, we show that non-linear wave equations which can be writ- 
ten as systems involving equations of the form □</> = (j>'V<j> and □</> = |V<^P 
are well-posed with scattering in (6 + 1) and higher dimensions if the Cauchy 
data are small in the scale invariant f.^ Besov space B""'^. This paper is the 
first in a series of works where we discuss the global regularity properties of 
general non-linear wave equations for all dimensions 4 ^ n. 



1. Introduction 

In this paper, our aim is to give a more or less complete description of the global 
regularity properties of generic homogeneous quadratic non-linear wave equations 
on (6 -f 1) and higher dimensional Minkowski space. The equations we will consider 
are all of the form: 

(1) = U{cl^,D4>) . 

Here □ :— —df + denotes the standard wave operator on R""*"^, and A/" is a 
smooth function of (j) and its first partial derivatives, which we denote by D(f). For 
all of the nonlincarities we study here, J\f will be assumed to be at least quadratic 
in nature, that is: 

Af{X,Y) = 0{\{X,Y)\^) , {X,Y) ^ 0. 

The homogeneity condition we require M to satisfy is that there exist a (vector) a 
such that: 

(2) AA(A'"0,A'"+iD0) = A'^+2AA(0, D0) , 

where we use multiindex notation for vector Af. The condition (0) implies that 
solutions to the system are invariant (again solutions) if one performs the scale 
transformations: 

(3) 0(-) - A'^0(A-) . 

The general class of equations which falls under this description contains virtually 
all massless non-linear field theories on Minkowski space, including the Yang Mills 
equations (YM), the wave-maps equations (WM), and the Maxwell-Dirac equa- 
tions (MD). We list the schematics for these systems respectively as: 



This work was conducted under NSF grant DMS-0100406. 



2 



JACOB STERBENZ 



(YM) 



UA = ADA + A^ , 



(WM) 



□ </« 



1^01 



2 



Uu = ADu 



(MD) 



UA = \Du\^ . 



The various values of a for these equations are (respectively) cr = 1, cr = 0, and 
(7 = 1). For a more complete introduction to these equations, see for instance 
the work [2] and For the purposes of this paper, will will only be concerned with 
the structure of these equations at the level of the generic schematics (YM)-(MD). 

The central problem we will be concerned with is that of giving a precise descrip- 
tion of the regularity assumptions needed in order to guarantee that the Cauchy 
problem for the system is globally well posed with scattering (GWPS). That is, 
given initial data: 



we wish to describe how much smoothness and decay (/, g) needs to possess in order 
for there to exist a unique global solution to the system with this given initial 
data. We also wish to show that the solutions we construct depend continuously 
on the initial data, and are asymptotic to solutions of the linear part of (Q. We 
will describe shortly in what sense we will require these notions to hold. 

Our main motivation here is to be able to prove global well-posedness for non- 
linear wave equations of the form in a context where the initial data may not 
be very smooth, and furthermore does not possess enough decay at space-like in- 
finity to be in L^. Also, we would like to understand how this can be done in 
situations where the equations being considered contain no special structure in the 
non-linearity. For instance, this is of interest in discussing the problem of small 
data global well-posedness for the Maxwell-Klein-Gordon and Yang-Mills equa- 
tions with the Lorentz gauge enforced instead of the more regular Coulomb gauge. 
This provides a significant point of departure from earlier works on the global 
existence theory of non-linear wave equations, which for the general case requires 
precise control on the initial data in certain weighted Sobolev spaces (see @]), or else 
requires the non-linearity to have some special algebraic or gauge structure which 
allows one to exploit some null form identities or apply an appropriate renormal- 
ization to the equation being considered (see ^3] and |12|). 

From the point of view of homogeneity, we are lead directly to considerations 
of the low regularity properties of the equations as follows: By a simple scal- 
ing argument^: one can see that the most efficient based regularity assumption 

^In conjunction with finite time blowup for large data. This phenomena is known to happen 
for higher dimensional equations with derivative non-linearities even in the presence of positive 
conserved quantities. 



(4) 



m 



9t0(O) = g , 
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possible on the initial data involves Sc — f ~ o" derivatives. Again, by scale in- 
variance and looking at unit frequency initial data^, one can see that if we are to 
impose only an smallness condition on the initial data which contains no phys- 
ical space weights, then Sc = ^ — cr is in fact the largest amount of derivatives we 
may work with. This leads us to consider the question of GWPS for initial data in 
the homogeneous Besov spaces B^^'P, for various values of p. In this work, we will 
concentrate solely on the case p — 1. This is the strongest scale and translation 
invariant control on the initial data possible, and will be crucial for the kind of 
non-linearities we work with here. In fact, it does not seem possible to push any 
type of global regularity for equations of the type |^ which contain derivatives in 
the non-linearity down to the scale invariant Sobolev space H'^" — B^"''^ unless 
the equations under consideration possess a great deal of special structure in the 
non-linearity. This has been done for the wave-maps equations (see |ll2j ^ and more 
recently for the Maxwell-Klein-Gordon equations in (6 + 1) and higher dimensions 
with the help of the Coulomb gauge (see [H]). 

In recent years, there has been much progress in our understanding of the low 
regularity local theory for general non-linear wave equations of the form . In the 
lower dimensional setting, i.e. when n = 2,3,4, it is known from counterexamples 
of Lindblad (see ^U]) that there is ill posedness for initial data in the Sobolev space 
ff**" , where sq ^ Sc + ■ Intimately connected with this phenomena is the failure 
of certain space-time estimates for the linear wave equation known as Strichartz 
estimates. Specifically, one docs not have anything close to an L^{L*) estimate in 
these dimensions. Such an estimate obviously plays a crucial role (via Duhamel's 
principle) in the quadratic theory. However, using the Strichartz estimates available 
in these dimensions along with Picard iteration in certain function spaces, one can 
show that the Lindblad counterexamples are sharp in there is local well-posedness 
for initial data in the spaces H'^ when Sc + < s (see for example [3]). 

In the higher dimensional setting, i.e. when the number of spatial dimensions is 
n = 5 or greater, one does have access to Strichartz estimates at the level of L^{L*) 
(see |3]), and it is possible to push the local theory down to where < e is 

arbitrary (see ^^). 

In all dimensions, the single most important factor which determines the local 
theory as well as the range of validity for Strichartz estimates is the existence of 
free waves which are highly concentrated along null directions in Minkowski space. 
These waves, known as Knapp counterexamples, resemble a single beam of light 
which remains coherent for a long period of time before dispersing. For a spe- 
cial class of non-linearities, known as "null structures", interactions between these 
coherent beams are effectively canceled, and one gains an improvement in the lo- 
cal theory of equations whose nonlinearities have this form (see for example [S], |S]). 

In both high and low dimensional settings, the analysis of certain null struc- 
tures, specifically non-linearities containing the Qo null from'^, has led to the proof 



That is initial data sets where the Fourier transform is supported in the unit frequency annulus 
I I < l?l< 2}. 

■^This is defined by the equation Qo{(f>,ip) = da(l>d°'ip. 
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that the wave-maps model equations* are well posed in the scale invariant £^ Besov 
space B^'^ (see |^ and [3). While the proof of this result is quite simple for high 
dimensions, it relies in an essential way on the structure of the Qo null form. In fact, 
there is no direct way to extend the proof of this result to include the less regular 
nonlinearities of the form or for that matter the Qij null forms^, which show 

up in the equations of gauge field theory. However, the high dimensional non-linear 
interaction of coherent waves is quite weak (e.g. giving the desired range of validity 
for Strichartz estimates), and one would expect that it is possible to prove local 
well posedness for quadratic equations with initial data in the scale invariant £^ 
Besov space without resorting to any additional structure in the nonlinearity. For 
n — 5 dimensions, it may be that this is not quite possible, although we provide 
no convincing evidence except for the fact that there is no obvious way to add over 
our localized estimates in that dimension in order to obtain a full set of estimates 
that works in all of space-time Fourier space. For n = 6 and higher dimensions, we 
will prove that in fact no null structure is needed for there to be well posedness in 
B^a,! 'jjjjg leads to the statement of our main result which is as follows: 



Theorem 1.1 (Global Well Posedness). Let 6 ^ n be the number of spatial dimen- 
sions. For any of the generic equations listed above: YM, WM, or MD, let (/, g) 
be a (possibly vector valued) initial data set. Let Sc = ^ — cr be the corresponding 
L^ scaling exponent. Then there exists constants < eo,C such that if 

(5) II if,9) lls^cixB^o-i,! «S eo , 

there exits a global solution ip which satisfies the continuity condition: 

(6) II V'llc(s=-i)nc(i)(s=c-i,i) C|| (/,.g) ||^a<„ixi3=c-i.i ■ 

The solution -tjj is unique in the following sense: There exists a sequence of smooth 
functions {fN,gN) such that: 

Jim II (/,5) - (/jv,5jv) IIb.o.ixs-c-i.i = 0. 

For this sequence of functions, there exists a sequence of unique smooth global so- 
lutions V'TV o/ CD with this initial data. Furthermore, the ipN converge to ip as 
follows: 

^lim^||V'-V'A'llc(B=c.i)nc(i)(s=<:-i.i) = 0- 

Also, '0 is the only solution which may be obtained as a limit (in the above sense) 
of solutions to ^ with regularizations of (/, g) as initial data. Finally, '0 retains 
any extra smoothness inherent in the initial data. That is, if (/, g) also has finite 
X H"~^ norm, for Sc < s, then so does ip at fixed time and one has the following 
estimate: 

(7) II V'llc(i/^)nc(i)(i/-i) «^ C\\{f,9)\\Hs^Hs-i ■ 



*Not the rough schematic we have hsted here, but rather equations of the form 

□</. = r(9i)Qo(</',0). 

^These are defined by Qij{<f),tp) = di<j),djip — dj4>diip. 
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In a straightforward way, the function spaces we iterate in aUow us to show the 
following scattering result without any extra work: 

Theorem 1.2. Using the same notation as above we have that there exists data 
sets {f^,g^), such that if is the solution to the homogeneous wave equation 
with the corresponding initial data, the following asymptotics hold: 

(8) ^lim ||V'+-V'llsac.ina«B-c-i.i = 0, 

(9) ^Um^ll V" - V'llsa=.ina«B»c-i,i = 0. 

Furthermore, the scattering operator retains any additional regularity inherent in 
the initial d.a.ta. That is, if {f,g) has finite norm, then so does {f^,g^), and 
the following asymptotics hold: 

(10) II - V'llij.natijo-i = 0' 

(11) ^hm^||V--V'll™,ij-i = 0- 



2. Preliminary Notation 

For quantities A and B, we denote hy A < B to mean that A ^ C ■ B for some 
large constant C. The constant C may change from line to line, but will always 
remain fixed for any given instance where this notation appears. Likewise we use 
the notation A ^ B to mean that ^ ■ B ^ A ^ C ■ B. We also use the notation 
A <^ B to mean that A ^ ^ ■ B foY some large constant C. This is the notation 
we will use throughout the paper to break down quantities into the standard cases: 
A ~ B, or A <C B, or i? <C A; and A < B, ov B <^ A, without ever discussing 
which constants we are using. 

For a given function of two variables {t, x) € TL x we write the spatial and 
space-time Fourier transform as: 

f{t,0= J e-2"«-/(t,x)dx, 

f{r,0 = j e-2-(-*+€-) f{t,x) dtdx . 

respectively. At times, we will also write 3^[f] = f. 

For a given set of functions of the spatial variable only, we denote by W{f, g) 
the solution of the homogeneous wave equation with Cauchy data (/, g). If F is a 
function on space-time, we will denote by W{F) the function W {F{0), dtF (0)). 

Let E denote any fundamental solution to the homogeneous wave equation, i.e., 
one has the formula UE = 5. We define the standard Cauchy parametrix for the 
wave equation by the formula: 

-F = E*F-W{E*F) . 
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Explicitly, one has the identity: 



For any function F which is supported away from the light cone in Fourier space, 
we shall use the following notation for division by the symbol of the wave equation: 

4f = E*F . 



Of course, the definition of i docs not depend on E so long as for F is supported 
away from the light cone; for us that will always be the case when we use this 
notation. Explicitly, one has the formula: 



47r2(T2 - |$|2) 



3. Multipliers and Function Spaces 

Let (/? be a smooth bump function (i.e. supported on the set \s\ ^ 2 such that 
= 1 for |s| < 1). In what follows, it will be a great convenience for us to assume 
that ip may change its exact form for two separate instances of the symbol ip (even if 
they occur on the same line) . In this way, we may assume without loss of generality 
that in addition to being smooth, we also have the idempotence identity = ip. 
We shall use this convention for all the cutoff functions we introduce in the sequel. 

For A e 2^, we denote the dyadic scaling of (p> by <^a(s) = fij)- The most basic 
Fourier localizations we shall use here arc with respect to the space-time variable 
and the distance from the cone. Accordingly, we form the Littlewood-Paley type 
cutoff functions: 

(13) SA(T,0=m(l(T,6l)-nA(K^'?)l) ' 

(14) cdT,0 = ^2d{\r\-\^\)-ViA\r\-\C\) • 

We now denote the corresponding Fourier multiplier operator via the formulas 
S\u = s\u and C^u = CdU respectively. We also use a multi-subscript notation 
to denote products of the above operators, e.g. Sx,d = SxC^ ■ We shall use the 
notation: 

(15) S\^»^d = Sx^s , 

to denote cutoff in an 0(d) neighborhood of the light cone in Fourier space. At 

times it will also be convenient to write Sx,d^, = Sx — Sx.,<d ■ We shall also use 
the notation S^^^ etc. to denote the multiplier Sx,d cutoff in the half space ±t > . 
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The other type of Fourier locaHzation which will be central to our analysis is the 
decomposition of the spatial variable into radially directed blocks of various sizes. 
To begin with, we denote the spatial frequency cutoff by: 

(16) PA(e) = ¥'2A(|^|)-^iA(lel) , 

with Px the corresponding operator. For a given parameter 6 ^ X, we now decom- 
pose P\ radially as follows. First decompose the the unit sphere 5""^ C M" 
into angular sectors of size x . . . x ^ with bounded overlap (independent of 
S). These angular sectors are then projected out to frequency A via rays through 
the origin. The result is a decomposition of supp{px} into radially directed blocks 
of size A X (5 X ... X (5 with bounded overlap. We enumerate these blocks and 
label the corresponding partition of unity by b'^g. It is clear that things may be 
arranged so that upon rotation onto the ^i-axis, each 6^ g satisfies the bound: 

(17) \d^bU ^ CjvA-^ , Idt'bU ^ CnS-^ . 

In particular, each B'^ g is given by convolution with an kernel. We shall also 
denote: 

Note that the operators S^d and Sx,<:d are only supported in the region where 
\r\ < 1^1 . 

We now use these multipliers to define the following dyadic norms, which will be 
the building blocks for the function spaces we will use here. 

(18) lluf 1 d^\\SxduE2 , ("classical" H'-^) 

(19) \\u\\y,=\-'\\DSxu\\li,^l-) , (Duhamel) 



(20) IhlU. =A^^ ^||5X<i«lliHi~) ■ (outer block) 

Notice that the (semi) norms ^ and Yx are only well defined modulo measures 
supported on the light cone in Fourier space. Because of this, it will be convenient 
for us to include an extra L°°(L^) norm in the definition of our function spaces. 
This represents the inclusion in the above norms of solutions to the wave equation 
with initial data. Adding everything together, we are led to define the following 
fixed frequency (semi) norms: 

(21) II u\\f, = (4,1 + Yx) n Sx {L^{L^)) ■ 

Unfortunately, the above norm is still not strong enough for us to be able to iterate 
equations of the form which contain derivatives. This is due to a very specific 
Low X High frequency interaction in quadratic non-linearities. Fortunately, this 
problem has been effectively handled by Tataru in jJS], based on ideas from [7] 
and What is necessary is to add some extra L^{L°°) norms on "outer block" 
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regions of Fourier space. This is the essence of the norm H2()|l above, which is a 
shght variant of that which appeared in |15j . This leads to our second main dyadic 
norm: 

(22) II u\\g, = [xl, + n) n Sx {L°^{L^)) n . 

FinaUy, the spaces we will iterate in are produced by adding the appropriate number 
of derivatives combined with the necessary Besov structures: 

(23) \\u\\l. = J^A^Ihll^^ , 

A 

(24) \\u\\g^ = E^II^IIg. • 

A 



Due to the need for precise microlocal decompositions, of crucial importance to 
us will be the boundedness of certain multipliers on the components l|18|) - H19|l of 
our function spaces as well as mixed Lebesgue spaces. We state these as follows: 



Lemma 3.1 (Multiplier boundedness). 

(1) The following multipliers are given by kernels: X^^VSx, S'<^ j^, S'^ 
and {Xd)S~^ S'^ ^ . In particular, all of these are bounded on every mixed 
Lebesgue space L'(L''). 

(2) The following multipliers are bounded on the spaces L'^{L^), for 1 ^ q ^ oo: 
S\,d and S\,,^d- 



Proof of Lemma \y.l\ (1). First, notice that after a rescaling, the symbol for the 
multiplier A^^VS'a is a C°° bump function with 0(1) support. Thus, its kernel is 
in with norm independent of A. 

For the remainder of the operators listed in (1) above, it suffices to work with 
(A(i)S'~^S'^jj. The boundedness of the others follows from a similar argument. We 
let denote the symbol of this operator cut off in the upper resp. lower half plane. 
After a rotation in the spatial domain, we may assume that the spatial projection 
of is directed along the positive axis. Now look at x'*'(s, 77) with coordinates: 

^1 = + 
^7' = e' . 

It is apparent that x^(*7 ^) has support in a box of dimension '--•Ax \p\d x . . . x 
\/Adx d with sides parallel to the coordinate axis and longest side in the r\\ direction 
and shortest side in the s direction. Furthermore, a direction calculation shows that 
one has the bounds: 
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Therefore, we have that yields an kernel. A similar argument works for the 
cutoff function x^, using the rotation: 

ri' = ■ 

□ 



Proof of Lemma \^~l\ (2). We will argue here for S\^d- The estimates for the others 
follow similarly. If we denote by K^{t^x) the convolution kernel associated with 
^\cV then a simple calculation shows that: 



where supp{ip} is contained in a box of dimension ^ A x . . . x A x d with sides 
along the coordinate axis and short side in the r direction. Furthermore, one has 
the estimate: 

This shows that we have the bound: 

WK^hliL-) < 1, 

independent of A and d. Thus, we get the desired bounds for the convolution 
kernels. □ 



As an immediate application of the above lemma, we show that the extra Z\ 
intersection in the G\ norm above only effects the ^ portion of things. 

Lemma 3.2 (Outer block estimate on Y\). For 5 < n, one has the following 
uniform inclusion: 

(25) Yx C Zx . 



proof of (|25|l . It is enough to show that: 

First, using a local Sobolev embedding, we see that: 
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Therefore, using the boundedness Lemma mi it suffices to note that by Minkowski's 
inequahty we can bound: 



< W S\,du\\L^L^-) 



□ 



The last hne of the above proof showed that it is possible to bound a square sum 
over an angular decomposition of a given function in L^{L^). It is also clear that 

this same procedure works for the spaces because one can use Minkowski's 

inequality for the £^ sum with respect to the cone variable d. This fact will be of 
great importance in what follows and we record it here as: 

Lemma 3.3 (Angular reconstruction of norms). Given a test function u and pa- 
rameter 5 ^ \, one can hound: 



(26) IEII^aVII^^.^^^^ 



< 



4. Structure of the F\ spaces 

The purpose of this section is to clarify some remarks of the previous section and 
write down two integral formulas for functions in the F\ space. This material is 
all more or less standard in the literature and we include it here primarily because 
the notation will be useful for our scattering result. Our first order of business is 
to write down a decomposition for functions in the -Fa space: 

Lemma 4.1 {F\ decomposition). For any u\ ^ F\, one can write: 

(27) u\ — Uy +Uyi/2+UY. , 

A,l 

where u-^^ is a solution to the homogeneous wave equation, U-^1/2 is the Fourier 
transform of an function, and uy^ satisfies: 

Furthermore, one has the norm bounds: 

(28) ^IIwaIIf^ =^ (^I|wxJIl~(l2) + ||u^i/2||^i +\\uy^\\y^^ ^ C|| ua ||f;, • 

We now show that the two inhomogeneous terms on the right hand side of H27|) 
can be written as integrals over solutions to the wave equation with data. This 
fact will be of crucial importance to us in the sequel. The first formula is simply a 
restatement of (I12II : 
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Lemma 4.2 (Duhamel's principle). Using the same notation as above, for any 
uyx, one can write: 

(29) uyjt) = - f \D,\-' sin {{t - s)\D,\) OuyAs) ds . 



Likewise, one can write the ■u„i/2 portion of the sum (I27II as an integral over 

modulated solutions to the wave equation be foliating Fourier space by forward 
and backward facing light-cones: 

Lemma 4.3 {X^^ Trace lemma). For any u^i/2, let u denote its restriction 
to the frequency half space < ±t. Then one can write: 

(30) u^,,^{t) = J e2"*(«±l^==l)u±^ ds , 

where ^ is the spatial Fourier transform of restricted to the s*'* translate of 
the forward or backward light-cone light cone in Fourier space, i. e. : 

^±A,s(0 = / Sir~sT\i\)^ir,0 dr . 
In particular, one has the formula: 

(31) / WutjL^ds < \\u%,\\i . 



5. Strichartz estimates 

Our inductive estimates will be based on a method of bilinear decompositions 
and local Strichartz estimates as in the work il^. We first state the standard 
Strichartz from which the local estimates follow. 

Lemma 5.1 (Homogeneous Strichartz estimates (see |3|)). Let 5 < n, a = ^^y^, 
and suppose u is a given function of the spatial variable only. Then */ ^ + 7 ^ § 
and ^ + 7 = f ^ 7; the following estimate holds: 

(32) lle^^'^^^l^-IP.^AiillLKLs) < X-'WP.^xuU^ . 

Combining the L {L "-3 ) endpoint of the above estimate with a local Sobolev in 
the spatial domain, we arrive at the following local version of H32|) : 

Lemma 5.2 (Local Strichartz estimate). Let 5 < n, then the following estimate 
holds: 

(33) ||g±2^»t|D.l5- ,u\\l2(i^o.. < X^d'^WB^ iu\\l2 . 
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Using the integral formulas lPt?)l and (jSOjlj we can transfer the above estimates to 
the Fa spaces: 

Lemma 5.3 (Fa Strichartz estimates). Let 5 < n and set a = ^^x'"- Then if 
I + ^ ^ ^ and ^ + 7 = § ~ 7; the following estimates hold: 

(34) \\Sxu\\L.(Ln ^ a^ii^iIf, , 



Proof of Lemma \5.!A It suffices to prove the estimate (j2U, as the estimate HH5|I 
follows from this and a local Sobolev embedding combined with the resumming 
formula H2t)|) . Using the decomposition (|4.1|) and the angular reconstruction formula 
(1^ . it is enough to prove for functions u 1/2 and uy_^. Using the integral 

formula ()30() . we see immediately that: 

I! u^^,. |U,(i.) IW e^'"*l^-l 4, ' 



< II li 1/2 II 1 



For the uy^ portion of things, we can chop the function up into a fixed number of 
space-time angular sectors using convolution kernels. Doing this and using 
to denote an operator from the set {/, i9i|-Da;|~^}, we estimate: 

II hiiL-} ^ II daUY^ Wl^L-) , 

a 

±,a •' 

<X''\-'Y1 J l|eT'^'^l^-li?on^/y,(s,a;)|U2 ds , 



^ XT' . 

□ 

A consequence of is that we have the embedding: 

XI, C i-(L2) . 

Using a simple approximation argument along with uniform convergence, we arrive 
at the following energy estimate for the F* and spaces: 
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Lemma 5.4 (Energy estimates). For space-time functions u, one has the following 
estimates: 

(36) II "llc(ff=)nc(i)(ij=-i) - 

(37) Il^^llc(s=)nc(i)(i3=-i) - II"IIg=- 

Also, by duality and the estimate (|36|l . we have that: 

(38) \S-'L\L') C XS-^X-^ C . 
This proves shows: 

Lemma 5.5 {LP' estimate for Y\). The following inclusion holds uniformly: 

(39) d^Sx^d{Yx) C ^2(^2) ^ 
in particular, by dyadic summing one has: 



6. Scattering 

It turns out that our scattering result, Theorem 11.21 is implicitly contained in 
the function spaces F'^ and ff*. That is, there is scattering in these spaces inde- 
pendently of any specific equation being considered. Therefore, to prove Theorem 
11.21 it will only be necessary to show that our solution to ^ belongs to these spaces. 

Using a simple approximation argument, it suffices to deal with things at fixed 
frequency. Because the estimates in Theorem II . 21 deal with more than one deriva- 
tive, we will show that: 



Lemma 6.1 {Fx scattering). For any function ux G there exists a set of initial 
data {f^,g^) G Px{L'^) 'X \Px{L'^) such that the following asymptotic holds: 

(40) ^Inn II ux{t) - W{ft,gt){t) H^ina^CL^) = » ' 

(41) lim \\ux{t)-W{f-,gl){t)\\^,^g^^^.,^ = 0. 



Proof of Lemma \6.1\ Using the notation of Section^ we may write: 

■"A = Uj^ + M+i /2 + ^^1/2 -I- UY^ , 

We now define the scattering data implicitly by the relations: 

/•oo 

W{f+,gl){t) = uj,^- \D^\-^ sin {\D^\{t ~ s)) DuyAs) ds , 
Jo 



W{fx,9xm = " 



.0 

/ \D,\~^ sin {\D,\{t-s)) any Js)ds 

J —OO 
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Using the fact that Ouy^ has finite L^{L'^) norm, it suffices to show that one has 
the limits: 

1™ II "^i/sW +""i/2W llffinOtfL^) = 0. 
Squaring this, we see that we must show the hmits: 



(42) hm \D,\u+ ,,{t) \D.Aul,,,{t) ^ 



(43) hm / dtU+,/,{t) dtu%,j,{t) = . 

We'U only deal here with the limit (|42|l . as the limit (|43|l follows from a virtually 
identical argument. Using the trace formula (|30|l along with the Plancherel theorem, 
we compute: 

{L.H.S.)^ 
By H31(l we have the bounds: 



ICP / 14.,+.. (01 • hL,(e)| ds.ds. 



< II 1 \\u^,,, II 1 

II V-1/2 II 2 " A ' 



This shows that the function: 



m) - ler J e™<,^+,,(c) <,,(o d^^irf52 , 

is bounded pointwise by an function uniformly in t. Therefore, by the dominated 
convergence theorem, it suffices to show that we in fact have that limt^±oo Ht = 0. 
To see this, notice that by the above bounds in conjunction with Fubini's theorem, 
we have that the integral: 



1^1' j<s,+JOuijO ds,, 

is m Ll^. The result now follows from the Riemann Lebesgue Lemma. Explicitly, 
one has that for almost every fixed ^, the following limit holds: 



^hm^i?t(e) = ^hrn^ \e J e'-''''ul^JO ^^i ^^^2 



□ 



7. Inductive Estimates I 

Our solution to will be produced through the usual procedure of Picard 
iteration. Because the initial data and our function spaces are both invariant with 
respect to the scaling any iteration procedure must effectively be global in time. 
Therefore, we shall have no need of an auxiliary time cutoff system as in the works 
Instead, we write directly as an integral equation: 

(44) <P = W^(/,5) + n-iW,Z?0) . 
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By the contraction mapping principle and the quadratic nature of the nonlinearity, 
to produce a solution to (|44|l which satisfies the regularity assumptions of our main 
theorem, it suffices to prove the following two sets of estimates: 

Theorem 7.1 (Solution of the division problem). Let 5 < n, then the F and G 
spaces solve the division problem for quadratic wave equations in the sense that for 
any of the model systems we have written above: YM, WM, or MD, one has the 
following estimates: 

(45) \\U-^M{u,Dv)\\g^. < ||7i||G=c||i;||G=c , 

(46) \\B-^Miu,Dv)y. < \\u\\G^4vy.+\\u\\F4v\\G^^ . 



The remainder of the paper is devoted to the proof of Theorem 17.11 In what 
follows, we will work exclusively with the equation: 

(47) = W{f,g) + D-\^V^) . 

In this case, we set Sc = The proof of Theorem 17. II for the other model equa- 
tions can be achieved through a straightforward adaptation of the estimates we 
give here. In fact, after the various derivatives and values of Sc are taken into ac- 
count, the proof in these cases follows verbatim from estimates ()49|l and (|5()|l below. 

Our first step is to take a Littlewood-Paley decomposition of □^^(mVu) with 
respect to space-time frequencies: 

(48) 0-\uVv) = ^□-i(5^,?/V5^,«) . 

fj.i 

We now follow the standard procedure of splitting the sum 148f) into three pieces 
depending on the cases /xi <C M2 ^ Mi: and ^2 Mi- Therefore, due to the £^ 
Besov structure in the F spaces, in order to prove both (|45|) and (|46|l . it suffices to 
show the two estimates: 



(49) \\0-\S,,uVS^,v)\\g, < X-'^il\\u\\Fj\v\\p^^^ ,Mi'-M2, 

(50) \\a-\S^uWS^v)\\G, < M'^lhllcjlf IIf. , m«A. 



Notice that after some weight trading, the estimates H45(l and (|46|l follow from H5()() 
in the case where 112 <C /ii. 

proof of (|^ . It is enough if we show the following two estimates: 

(51) \\Sx{S^,,u\/Sf,^v)\\Li{L--) < I^I\\u\\f^,\\v\\f,,^ ,Aii-Ai2, 

(52) \\Sxa-\Sf,,uVSf,^v)\\L^i^L2) < A"Vl"l|u||F,,Jh||F^, ,/ii~M2- 
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In fact, it suffices to prove H51fl . To see this, notice that one has the formula: 

[Sx,D~^]G = W{E*SxG)-SxW{E*G) . 
Thus, after multiplying by 5*;^, we sec that: 

5A[5A,n"']G = Px{W{E*SxG)-SxW{E*G)) , 
= WiE * SxPxG) - SxW{E * PxG) , 
= Sx[Sx,a-^]PxG . 

Therefore, by the (approximate) idenipotence of Sx one has: 

SxO-^G = SaD-^^aG + S'a [5A,n"'] G , 

= Sxa-'SxG + Sx[Sx,a-^]PxG . 

Thus, by the boundedness of 5a on the spaces L°°{L'^) the energy estimate, one 
can bound: 

||5An-iG|Uoo(i.) < A-i(|15AG|Ui(i2) + |lPAG|Ui(i2) + |15APAG|Ui(i2)) . 



We now use the fact that the multipliers Sx and Pa are both bounded on the space 
L^(L'^) to reduce things to the estimate: 

II S'^iM VS'p^i; ||li(l2) < II S'piU ||l2(l4)|| VS'^^-y ||l2(l4) , 

Tl-2 Tl + 2 

II "IIf^JI w||f^, ■ 

Taking into account the the bound ^1 ~ the claim now follows. □ 

Next, we'll deal with the estimate (|5U|I . For the remainder of the paper we shall 
fix both A and fi and assume they such that /i <C A for a fixed constant. We now 
decompose the product Sx{Sf^uV Sxv) into a sum of three pieces: 

Sx{Sf,uWSxv) = A + B + G , 

where 

A = S'a(S'^u VS'a.c^^.u) , 

B = Sx,ctisi»{Sf^uVSx,,<ctJ.v) , 

C = Sx,,<ct^{Sfj,uVSx,,<ctJ.v) ■ 

Here c is a suitably small constant which will be chosen later. It will be needed to 
make explicit a dependency between some of the constants which arise in a specific 
frequency localization in the sequel. We now work to recover the estimate (15011 for 
each of the three above terms separately. 
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proof of (|5()|l for the term A. Following the remarks at the beginning of the proof 
of l^nj, it suffices to compute: 

II S'a(S'^u VS'a,c^^,i;) < A II S"^?/ 11^2(2,^)11 S'a,c^<5,u II L2(i2) , 

< A/xV||u||i.^^ M^^lhllF, , 

< c-iAAiT^||u||Fj|«||i=^;. . 

For a fixed c, we obtain the desired result. □ 



We now move on to showing the inclusion H5U|I for the B term above. In this 
range, we are forced to work outside the context of L^{L^) estimates. This is the 

reason we have included the L'^{L'^) based X^-,^ spaces. This also means that we 
will need to recover Z\ norms by hand (because they are only covered by the Y\ 
spaces). However, because this last task will require a somewhat finer analysis than 
what we will do in this section, we contend ourselves here with showing: 



proof of the ^ n Sx{L°°{L'^)) estimates for the term B. Our first task will be deal 
with the energy estimate which we write as: 

For G supported away from the light-cone in Fourier space, we have the identity: 
SxO-'^SxG = S-^SxG-W{S-^PxSxG) . 

Therefore, by using the energy estimate for the X^ ^ space, this allows us to esti- 
mate: 

WSxa-^SxGU^^L-} < \\S'^SxG\\l^^l-) + \\W{S-^PxSxG)\\l^^l2^ , 
< II ^^^SxG ||l=°(l2) , 
h( „ , 



< \\S-'SxG\ 



Therefore, we are left with estimating the term B in the X^ ^ space. For a fixed 
distance d from the cone, we compute that: 

\\ ^~^Sx,dSx,ciJ.!i»iSfj.uySx,m<ctJ.v)\\L^(L^) < d^^ || ||l2(L~) || 'S'aW ||Loo(i2) , 

< d-V^I|f||Fjk|k. . 

Summing d^ times this last expression over all cfi ^ d yields: 
d^W Sx,dSx,ctisi,{Sf^uV Sx^,<ctiv) ||l2(l2) 



fi\2 S—1 

M " \\v\\fA\u\\f^ 



c^^d 

For a fixed c we obtain the desired result. 

□ 



18 



JACOB STERBENZ 



8. Interlude: Some bilinear decompositions 

To proceed further, it will be necessary for us to take a closer look at the ex- 
pression: 

(53) Sla{Si,uVSx,.<c^v) , cfi^d, 

as well as the sum: 

C = Sx,.<c^^{S^,uWSx,.<c^^v) = Ci + Cu + Cm , 



where 



d<cfj, 

Cii = ^ Sx,.<d{Sn,.^dU VSx,dv) , 

• <mm{c/Lt,(i} 

{S1J.4U VSx, • <min{ciJ,4}'^ ) 



We'll begin with a decomposition of C/ and Cu. The C/// term is basically the 
same but requires a slightly more delicate analysis. All of the decompositions we 
compute here will be for a fixed d. The full decomposition will then be given by 
summing over the relevant values of d. Because our decompositions will be with 
respect to Fourier supports, it suffices to look at the convolution product of the 
corresponding cutoff functions in Fourier space. In what follows, we'll only deal 
with the Ci term. It will become apparent that the same idea works for C//. 
Therefore, without loss of generality, we shall decompose the product: 

To do this, we use the standard device of restricting the angle of interaction in 
the above product. It will be crucial for us to be able to make these restrictions 
based only on the spatial Fourier variables, because we will need to reconstruct 
our decompositions through square-summing. For (t',^') e swpp{s^,^^} and 
(r, ^) e supp{sj,^^} we compute that: 



o{d) = \y + T\~\e + ^\\ , 

= |±|^'| + |^|+0(d)|-|^' + ^|| , 
= 0{d) + \±\^'\ + \^\\-\^' + ^\ . 



Using now the fact that d < cfj, and fj, < cX to conclude that |^'| 
1^1 A , we see that one has the angular restriction: 



fj, and 



±ri + 1^1-1^ + ^1 = 0{d). 
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In particular we have that 0±^'.^ < ^Jj^ . This aUows us to decompose the 
product (|54|l into a sum over angular regions with 0(^/ — ) spread. The result is: 



Lemma 8.1 (Wide angle decomposition). In the ranges stated for the Ci term 

above, one can write: 

(55) 

\bJl^uJ2\~(d/ 

|ui-W3|~(d//i)^ 

for the convolution of the associated cutoff functions in Fourier space. 



We note here that the key feature in the decomposition H55|l is that the sum is (es- 
sentially) diagonal in all three angles which appear there (wi, ^2, wa). It is useful 
here to keep in mind the following diagram: 



Figure 1. Spatial supports in the wide angle decomposition. 



We now focus our attention on decomposing the convolution: 

^^^^ ^A,«$Cmm{c/i,d} ^"^/^id * '^A,»5Cm'in{c/^,d} ^ 

If it is the case that d ^ pi , then the same calculation which was used to produce 
(|55|l works and we end up with the same type of sum. However, if we are in the 
case where d ~ /i , we need to compute things a bit more carefully in order to 
ensure that we may still decompose the multiplier s^ ^ using only restrictions in 
the spatial variable. To do this, we will now assume that things are set up so that 
cfj, <^ d . It is clear that all the previous decompositions can be made so that we 
can reduce things to this consideration. If we now take {t',^') G ^'^J-PpI^^ j} 
and (t, £_) e SMpp{s^,^„j„{^^^}} , we can use the facts that r' = 0^{d) ±\£_'\ , 
T = 0{cfi) + 1^1 , and |^| ^ fi to compute that: 
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OM= |T' + r|-|r + el 
(57) = OT(d)±|e'|+OM + |e|-|e' + ei 



where the term 0^{d) in the above expression is such that |OT(d)| d . In fact, 
one can see that the equahty (|^ forces ±0'^{d) < on account of the fact 
that ±(±|^'| + |^|-|^' + ^|) >0 and the assumption \0{c^i) + 0^{d)\ ^ d . In 
particular, this means that we can multiply ^ in the product (|56|l by the cutoff 
S|t|<|^I without effecting things. This in turn shows that we may decompose the 
product (|56|l based solely on restriction of the spatial Fourier variables, just as we 
did to get the sum in Lemma FS. II 

We now return to the C/ term. For the sequel, we will need to know what the 
contribution of the factor S'A,»^di' to the following localized product is: 

Using Lemma l8. II we see that we may write: 



where \uji — ±uj2\ ~ \'-^3 — i'^2| ~ • However, this can be refined significantly. 
To see this, assume that the spatial support of s^"^ lies along the positive axis. 
We'll label this block by b'^^ i . Because we are in the range where y/JId ^ ^/Xd , 

we see that since for any ^ G supp{b^^ d i) ^'^'^ ^' ^ ^"PH'{*^\<(i} ^'^'^ 
^ + ^' must belong to supp{b^^ i } , we in fact have that ^ itself must belong 

A, (Ad) 2 

to a block of size A x ^fXd x . . . x \/Xd . This allows us to write: 



Lemma 8.2 (Small angle decomposition). In the ranges stated for the Cj term 
above, we can write: 

(59) s^^d(s^,.<;d * SA,.^d) = s^'d(<'.s:d * ^^f.sSd) • 

where \u!i ~ cu^l ^ , and \u!i — ±uj2\ ^ ■ 

It is important to note here that if one were to sum the expression H59(l over oji, 
the resulting sum would be (essentially) diagonal in wa, but there would be many 
LOi which would contribute to a single lu2. This means that the resulting would not 
be diagonal in ui2 as was the case for the sum (|55|l . It is helpful to visualize things 
through the following figure: 



Our final task here is to mention an analog of Lemma HOI for the term (|53ll . Here 
we can frequency localize the factor 5'A,.<cAt in the product using the fact that one 
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/id 




Figure 2. Spatial supports in the small angular decomposition. 



has fJ, c 2 y/Xd . The result is: 



Lemma 8.3 (Small angle decomposition for the term B). In the ranges stated for 
the B term above, we can write: 



(60) 



A, (Ad)' 



where |wi — W3I 



Finally, we note here the important fact that in the decomposition H60() above, the 
range of interaction in the product forces d ^ n . This completes our list of 
bilinear decompositions. 



9. Inductive Estimates II: Remainder of the Low x High^ High 

FREQUENCY INTERACTION 

It remains for us is to bound the term B from line (|53|l in the Z\ space, as well 
as show the inclusion H50|l for the terms C/ - Cm from line H54() . We do this now, 
proceeding in reverse order. 



proof of estimate (|5U|I for the Cm term. To begin with we fix d. Using the re- 
marks at the beginning of the proof of we see that it is enough to show that: 

(61) II ^A. 

^ A(^||5-,u||i,(i^) ) \\v\\f, . 
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To accomplish this, we first use the wide angle decomposition, H55|l . on the left 
hand side of H61(l . This allows us to compute, using a Cauchy-Schwartz, that: 

II ^A, d}{S,^,dU VS'A,,<™„{c^,d}'y) ||Li(i=) ' 

~ X! II IliHi"") ■ II ^-^^^f da'^A,.<mm{c/i,d}^') IIl~(L2) , 

|w3±aJ2|~(<i/M)^ 

^ ^(Ell^AVIlin^-)) Ml^'^ ■ 

Summing over d now yields the desired estimate. □ 



proof of (|50|l for the C// term. Again, fixing d, and using the angular decomposi- 
tion lemma IHTTl we compute that: 

II »S'A,»<ci(>5'^t,» ^du WSx^dv) ||li(l2) , 

E ii'5Xk<i"iu^(i-)-ii^rAf^)^^'''^^"^^(^^) ' 

|a;3±Ci;2|~(c?//j,) 2 

< A (^E II ^X.^d" lli^(L~)) II ^A.d^; WmL^) , 
^ V 2 ( - j ||w||fJ|w||f;. . 

This last expression can now be summed over d, using the condition d < cfj,, to 
obtain the desired result. □ 



proof of (|5U|I for the Cj term. This is the other instance where we will have to rely 
on the ^ space. Following the same reasoning used previously, we first bound: 

II ^^^S\^d{Sf_,^,i:dU S/S\^,^dv) ||l2(^2) , 

^ d-' E \\s;:,^d^\\LHL^)-\\B-^ ,s,,,^dv)h^^L2), 

|aJ3±CJ2|~(d//J.)2 

< d''(T.\\s:i,.<du\\hiL^i] •(Eii^rw^)i'^A..^'^^)iii-(i^)) ' 



1 n-2 



< d 2^ 2 ( _ ) ||w||fJ|«||f^ 



Multiplying this last expression by di and then using the condition d < cfi to sum 

- 

over d yields the desired result for the ^ space part of estimate l|3n|) . It remains 
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to prove the Zx estimate. Here we use the second angular decomposition lemma 
to compute that for fixed d: 



sidV) ||l1(L~) 



\cJl±(^2~(<i/A')' 



Tl-5 71-5 

d\~ fd" 



< 



^-2 71-2 , 



2 A 2 ||u||i.J|w|| 



Multiplying this last expression by A 2 and summing over d using the condition 
d < X,iJ yields the desired result. □ 



proof of the Z\ embedding for the B term. The pattern here follows that of the last 
few lines of the previous proof. Fixing d, we use the decomposition Lemma 18.31 to 
compute that: 



< 



ll,i(L° 



< d ^||S'^u||L2(io=) • ( ^||B^^^^^iS'A,.s:cpu||^2(i^) 



< 



& 



5 /d\~ 
A 



m"^ a"^ ||u||fJ|w|| 



Multiplying the last line above by a factor of A 2 and using the conditions d < X 
and c/j, < d < 11 , we may sum over d to yield the desired result. □ 
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